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Abstract 

In this paper we obtain general conditions under which the wave 
equation is well-posed in spacetimes with metrics of Lipschitz regu¬ 
larity. In particular, the results can be applied to spacetimes where 
there is a loss of regularity on a hypersurface such as shell-crossing 
singularities, thin shells of matter and surface layers. This provides a 
framework for regarding gravitational singularities not as obstructions 
to the world lines of point-particles, but rather as obstruction to the 
dynamics of test fields. 


1 Introduction 

An important requirement of any classical physical theory is that given suit¬ 
able initial data one can determine the evolution of the system. Within the 
theory of general relativity the concept of global hyperbolicity [T] therefore 
plays a key role. Mathematically, a spacetime region Af is said to glob¬ 
ally hyperbolic if the causality condition is satisfied, and for any two points 
p,q £ Af the causal diamond J~^{p) n J~{q) is compact and contained in 
N [2]- However from the physical point of view the important property of 
smooth globally hyperbolic spacetimes is that the evolution of the Einstein 
equations is well defined (see e.g. [3] for details). 

One way of interpreting the compactness of J~^{p) 0 J~{q) is that this set 
“does not contain any points on the edge of spacetime, i.e. at infinity or at a 
singularity” O §6.6]. In this context the regularity usually considered in the 
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region M is that the metric must be (7^’^(also denoted This criteria is 

based on the existence and uniqueness of geodesics and the characterisation 
of gravitational singularities in terms of geodesic incompleteness. Moreover 
this regularity is the threshold where causal theory for smooth metrics and 
rough metrics agree mis]- 

However, in many realistic astrophysical situations one would like to solve 
Einstein’s Field Equations with lower regularity. For example, when one 
models a jump in density at the boundary of a star. This led Clarke [B] 
to define a notion of generalised hyperbolicity directly in terms of the local 
existence and uniqueness of the wave equation. 

Clarke’s definition of generalised hyperbolicity was motivated by two things. 
Firstly there are a number of spacetimes in which points are removed due to 
the presence of weak singularities and are therefore not globally hyperbolic 
but are still of physical interest. These include spacetimes with thin shells of 
matter [7j , impulsive gravitational waves [S] and shell-crossing singularities 
[9] . The second motivation was that of using test fields (given by solutions 
of the wave equation) rather than test particles (given by solutions of the 
geodesic equation) to probe the structure of singularities. Some work to 
examine the physical effect of gravitational singularities has been done by 
using a 3-parameter family of test particles (see for example m), however 
an advantage of using test fields is that the behaviour of the naturally de¬ 
fined energy-momentum tensor of the field gives a direct measurement of 
the physical effect of the singularity which is not so easy when considering 
families of test particles. 

There have been previous approaches which study the nature of the singu¬ 
larities using test fields (see e.g. [m [121 [13 E] ). In all these approaches 
one considers self-adjoint extensions of the (spatial) Laplace-Beltrami oper¬ 
ator and applies appropriate boundary conditions at the singularity. All this 
work has focused on the case of static spacetimes and furthermore most of 
these self-adjoint extensions are in (except [13] who consider finite energy 
solutions). 

The approach in this paper is to consider singularities as interior points of 
a spacetime with low regularity rather than regarding the singularity as a 
boundary and applying boundary conditions. A natural condition in this 
context is to require the solutions to lie in the Sobolev space which 
ensures that the energy momentum tensor is well-defined as a distribution. 
In this paper we will therefore look at conditions on the spacetime and initial 
data which give solutions with this regularity. 

In particular, we will look at the co-dimension one case which is relevant 
for example to junction conditions with jumps in the density, impulsive 

function / on an open set U of R" is said to be Lipschitz or (7°’^ if there is some 
constant K such that for each pair of points p,q £U, |/(p) — fiq)\ ^ K\p — q\, where |p| 
denotes the usual Euclidean distance. We denote by those functions where the fcth 
derivative is a Lipschitz function. 
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gravitational waves and brane-world cosmologies. In a previous paper we 
looked at the case of co-dimension two singularities such as cosmic strings 

US]. 

In previous work Grant et al m showed the existence of generalised so¬ 
lutions to the wave equation on singular spacetimes. This involved replac¬ 
ing the singular metric by a 1-parameter family of smooth metrics. The 
1-parameter family of solutions of the corresponding wave equations then 
describes a generalised solution. Using this approach it was possible to 
demonstrate the existence of a generalised solution to the wave equation for 
a wide class of metrics which included metrics with components bounded al¬ 
most everywhere. A downside of this approach is that it is not easy to relate 
the generalised solution to a classical weak solution of the equation. In our 
approach we will utilise the vanishing viscosity method (see e.g. Evans HD) 
to show the existence of a weak solution of the wave equation. This involves 
approximating the wave equation by a system of second order parabolic 
equations with a parameter e corresponding to the viscosity. Like Grant 
et al we again obtain a 1-parameter family of solutions but we are able to 
utilise results from parabolic regularity theory to have better control over the 
solutions and their convergence. This enables us to show the 1-parameter 
family converges to a weak solution of the zero viscosity equation which 
corresponds to an solution of the wave equation. The basic method we 
use follows that of Evans m but the details differ and our results are also 
distinct from from his since we assume less regularity and as a result our 
solutions also exhibit less regularity. 

The plan of the paper is as follows. In §2 we introduce the general setting for 
the problem and state the main theorems. We do not impose the ultrastatic 
condition used in [18] and write the wave equation as a first order system 
(see e.g. m)- This enable us to work with the energy of the first order 
system which corresponds to an energy of the second order system. In 
this setting we obtain existence of the solutions in time dependent scenarios 
without the need for any symmetries. §3 contains the proofs of the theorems 
while §4 gives examples of how these can be applied to various examples such 
as a discontinuity across a hypersurface, impulsive gravitational waves and 
brane-world cosmologies. 

2 The main results 

The interest in co-dimension one singular submanifolds covers a variety of 
different interesting physical phenomena such as, surface layers [7], impul¬ 
sive gravitational waves |H], and shell-crossing singularities [2U| all of which 
fall outside the class of smooth globally hyperbolic spacetimes. Moreover, 
the mathematical analysis by Geroch and Traschen m of what now is called 
the class of Geroch-Traschen metrics and the subsequent analysis by Stein- 
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bauer and Vickers m using generalised functions gives co-dimension one 
singular submanifolds a robust mathematical background. In addition, re¬ 
cent proposals in semi-classical gravity and quantum gravity m suggest 
that the metric near the event horizon may present some loss of regularity. 
In this section we present techniques to prove local well-posedness of the 
wave equation in spacetimes with co-dimension one singularities subject to 
certain condition on the metric. 

Clarke used the term generalised hyperbolicity to describe the situation where 
one has a unique solution to the wave equation. However we want to impose 
the slightly stronger condition of well-posedness and in addition want to 
emphasise the role of the wave equation. See |24j for a discussion of this and 
related terms. 

2.1 The general setting 

The geometric setting is a region S[o,t] = [OiT] x S, where S is either a 
compact closed n-dimensional manifold or an open, bounded set of M” with 
smooth boundary dTi. In what follows, for simplicity, we will only consider 
the former case. The proof in the latter case follows by replacing by 

Hq(T,) and using the volume form given by dx^. 

Rather than considering the particular case of a spacetime with a singular 
hypersurface, where the regularity of the metric drops below we will 
consider a rough spacetime, where the spacetime metric gab is only Lips- 
chitz. We will show that in this situation one has well-posedness (in a sense 
made precise below) of the wave equation with weak solutions of regularity 
2 ’]). In order to do this, we will reformulate the wave equation as 
a first order symmetric hyperbolic system and look for solutions of this 
system. 

We therefore start by considering the first order initial value problem 

Lu = A^dtu +A^diu +Bu = F (1) 

u(0,-) = uo(-) (2) 

In the above we have employed the Einstein summation convention where i, 
j, k, etc. range over 1, 2,... , n. The unknown u and the source term F are 
both valued functions on while and B are N x N matrix 

valued functions on We will assume that A^ and Ad are symmetric 

and that in addition AA is positive definite. 

In order for such a system to correspond to the wave equation given by a Lip- 
schitz metric, we will require that A^ and Ad have bounded first derivatives 
and that B is bounded, so that we require 
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In the analysis below we will be working in spaces such as L^(S). Rather 
than defining this in terms of a particular coordinate system on S, we will 
introduce a background Riemannian metric hij on S and let Vh be the corre¬ 
sponding volume form. We then define to be the space of real valued 

functions ^ on S such that < oo and we denote the associated in¬ 
ner product by if,g)L^CE) = Note that since S is compact Uh is 

bounded from below and above. Furthermore if S is parallelizable (which 
for simplicity we will assume) there is no loss of generality in taking hij 
to be the flat metric. Note that in the three dimensional case, which is 
most relevant to applications in general relativity, it is enough for S to be 
orientable for it to be parallelizable. 

For the case of vector valued functions v on S, we dehne L^(S,M'^) to 
be those v such that v-viy^ < oo. The corresponding inner product on 
L^(E,R^) is then given by 


(v,w)i2(s,RJV) = 

Where there is no risk of confusion, we will write both the inner product 
in LP‘{Ti) and in simply as (•, ■)i ^2 The Sobolev spaces 

etc. are dehned in a similar manner ( see[T7] §5.2, |19|f . 

We also make use of the function space L^(S[o,r]) which is defined by re¬ 
quiring that functions are square integrable on [0, T] x S with respect to the 
volume form dt However in the analysis below it is often convenient to 
think of a function v(f, x) as a map from [0, T] to a function v(t)(-) of x € S 
given by v(f)(x) = v(f,x). For example L^(0, T; L^(S, is the space of 
functions 


v:[0,r] ^ (3) 

t ^ v(f) (4) 

such that v(f) G and 

^ (v,v)L 2 (s,MJV)dt < oo (5) 

When thinking of v in this way, we will denote the time derivative by v. 
We also dehne the following spaces which will be used in §3.4. 

W^’P(S[o,t],K^) := {w G W^’P(S[o,t],K^) : such that w(r, •) = 0} 

In particular, we will be interested in the following case 

:= {w G : such that w(r, •) = 0}. 
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2.2 Weak solutions and the main theorems 

We will be looking for weak solutions of the initial value problem ([1]). To 
motivate the definition we proceed as follows. Given a standard solution 
u of the initial value problem, we may first take the dot product of equation 
m with a smooth M^-valued function v with support in [0, T) x S and then 
integrate over x and t to obtain 

(Tu, v)£ 2 ('s,]Rjv)dt = ^ (F, (6) 

Integrating the left hand side by parts with respect to x and t we obtain 

^ (F,v)i2(s,RiV)(it (7) 

where L* is the formal adjoint of L defined below and the second term on 
the right hand side (RHS) comes from the t = 0 boundary term when we 
integrate by parts with respect to t. This approach results in the following 
definition: 

Definition 1 (^Weak Solution) We say a function: 

u e L2(0,r;L2(S,M^)) 

is a local weak solution of the initial value problem m provided that: For 
all V G with supp(v) C [0,T) x S 

(u, L*v) 2 , 2 (s,]RJV)dt = ^ (F, v)i 2 (s,MAf)dt + (^°(0)uo, . 

’ ( 8 ) 

This definition of a weak solution is the classical one used by Friedrichs [25] 
but differs from the one used by Evans m, who does not integrate with 
respect to t. Note also that the formal adjoint is defined with respect to 
Uh. So in the case where we use a general Riemannian metric, there are 
additional terms involving the derivatives of in the expression for L* 
compared to the flat case. The explicit expression is: 

L*v := -5t(24°v) - di{Wv) + - tIaW (9) 

where T*-^ are the connection coefficients of the smooth Riemannian metric 
hij. 

In order to prove uniqueness and well-posedness of the initial value problem, 
we will need to control the Lf size of the solution. This motivates the 
following definition. 
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Definition 2 (Tlegular Weak Solutionj We say a weak solution 
u G L^(0, T; M^)) is regular if u satisfies the energy estimate 

ll^llL2(o,r;L2(s,lR^)) ^ C" ^||uo||i2(s,MJV) + 1') I li2('s,]RJV)dt^ (10) 

We may now state our main result concerning solutions of low-regularity 
symmetric hyperbolic systems. 

Theorem 1 Given the linear symmetric hyperbolic system: 

Lv = A^dtu +A^diu +Bu = F (11) 

u(0,-) = uo(-) (12) 

where A*, i? and F are as above, and the initial data uq is in T^(S, M^). 
Then there exists a unique regular weak solution u G L^(0,T;L^(S,M^)). 
Furthermore this solution is stable in the sense that the solution depends 
continuously on the norm of the initial data in L^(S,M^) and the norm of 
the source function in L^(S[o^ 7 ’],M'^). 

We may now use the above result to establish the following theorem for the 
wave equation. 

Theorem 2 Let gab:9°'^ be in and f in L^(S[o^t])- Given initial data 
{uo,ui) G X LfiiYi) then the system 


\AgU + mfiu 

= f 

(13) 

m(0, •) 

= Uo 

(14) 

dtufi), ■) 

= Ul 

(15) 


has a unique stable solution u G Moreover, the corresponding 

energy-momentum tensor Tab[u] is in L]^^(S[o,r])- 

Note that the above result is similar to the one obtained for the homogeneous 
wave equation in [2B]. However, we have extended the results to the more 
general case that includes mixed space and time derivatives,included a source 
function and provided a different proof. 

We would like to mention that one can not go further without paying a price. 
As shown by Colombini et al. m there are examples of wave equations with 
coefficients depending only on time of Holder regularity!! with exponent 
a < 1 with no distributional solution. However, they proved well-posedness 
by moving from Sobolev spaces to Gevrey spaces. These results have been 

^ A function / on an open set U of R" is said to be Holder or (7°’“ if there is some non 
negative constant K such that for each pair of points p,q GlA, |/(p) — f{q)\ ^ K\p — q\°‘, 
where \p\ denotes the usual Euclidean distance. 
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further extended to coefficients that depend smoothly in the space variable 
but are log-Lipschitz (LL) H regular in time [5H]. In both [27] and [25| the 
structure of the second order part of the operator considered has the special 
form 



In |2S] mixed terms in time and space were allowed and the regularity of 
the coefficients was LL in space and time. However, the local well posedness 
results obtained are in Sobolev spaces with |s| < 1 and therefore the 
energy momentum tensor of the solutions is not integrable. Finally, one 
can explore wave-type equations with very rough coefficients such as in m 
where one is lead to the necessity of weakening the notion of a solution to 
the Cauchy problem and enlarging the allowed class of solutions. 

3 Proof of the main theorem 

3.1 Outline of the proof 

The proof of Theorem |2] uses the vanishing viscosity method described in 
§7.3 of Evans m- Note however that Evans assumes that the and A® 
have greater regularity than we do and as a result is able to obtain a solution 
with greater regularity. This explains why our definition of a weak solution 
has to differ from his. However the essence of the proof is essentially the 
same. It consists of the following steps: 

1. First, we approximate the problem © by the system of parabolic 
initial value-problem on E[o,r] given by 

dtu^-eAhU^ + {A\^A^diU^ + {A^)-^Bu^ = (7l°)-^F (17) 

u^(0,x) = p"*(uo(x)) 

where {(p'^)} G (0,1] is a family of mollifiers. Here Ah is the Laplace- 
Beltrami operator on T, associated with the smooth background Rie- 
mannian metric hij. By adding in the second order Laplace-Beltrami 
terms we obtain a system with smooth principal symbol. We may 
then use classical methods of parabolic regularity theory to obtain a 
solution with better analytic properties than the original hyperbolic 
system. 

function / on an open set U of R" is said to be log-Lipschitz if there is some 
non negative constant K such that for each pair of points p,q & ht, |/(p) ~ /(9)l ^ 
K\p — q|| In \p — g||, where |p| denotes the usual Euclidean distance. 




2. Second, we obtain the following uniform energy estimate 

^||uo||j;,2(s,RJV) + ^ I |F(t,') I lL2(s,R^)C?i 

(18) 

where C is independent of e. 

3. Third, we take the limit e —)■ 0 and show convergence in an appropriate 
weak sense to a regular weak solution as dehned above. 

4. Fourth, using the energy inequality (|18p we show uniqueness and sta¬ 
bility. This concludes the proof of Theorem [T] 

5. Fifth, we rewrite the wave equation as a symmetric hyperbolic problem 
and show that for a Lipschitz metric the corresponding L satisfies the 
conditions of Theorem [TJ 

6. Sixth, we show that the solution of the wave equation obtained via 

the symmetric hyperbolic problem is in This concludes 

the proof of Theorem [21 

3.2 Approximate solutions and energy estimate 

The results we obtain make extensive use of the vanishing viscosity method. 
As explained above, the first step is to show that there exist suitable solu¬ 
tions to (|17l) . This step follows directly from the work of Evans ( |17j , Th. 1 
§7.3). 

Proposition 1 ( Existence of Approximate solutions) For each e > 0, 
there exists a unique solution of f77lj with G and 

G L2(0,r;T2(S,M^)). 

Proof. This is a variant of a standard result for parabolic systems. Following 
Evans set X = L°°(0, T; F7^(S, M^)) and then for each w G A, consider the 
linear system 

dtu^-eAhU^ = (A°)-^5w +(A°)-^F (19) 

u'^(0, x) = Uq(x) (20) 

where Uq(x) = p^{x) * (uo(x)). Notice that the system is formed of N scalar 
parabolic equations of the form dtv — eA^v = f. The coefficients are now 
all smooth and the only loss of regularity comes from the source term on 
the RHS of (fTU|) . However as this is bounded in L^(0, T ; T^(S, M-^)), we can 
apply standard results (see e.g. m Th. 5 §7.1) to show that there is a 
unique solution G L^(0, T; F7^(S, M^)) with u'^ G T^(0, T; L^(S, M^)). 

In the same manner we can choose w G A and find u*^ that solves 
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- eAftU^ 
u^(0, x) 


( 21 ) 

( 22 ) 


= -(A°)-^^*9iW-(A°)-^5w +(A°)-^F 

= u3(a;) 

Subtracting u*^ — u^, we find u*^ = u*^ — solves 

atu"-eA;,u" = -{A^)-^A^diW-{A^y^Bw (23) 

u"(0,x) = 0 (24) 

where w = w—w. Using standard energy estimates for solutions of parabolic 
equations we have that u satishes: 




sup ||u (i)|Rjv) 

OsgtsJT ^ ' 

C{T, e) (||(A0)-i(^*9iW + 

< C{T,e) ^^sup^||w(t)||^i(j._jg;v)^ 


Thus 


u 


lL°°( 0 ,r;(E,R^)) ^ C'(r, e)||w||2,cxj(o^7’;(s 


Therefore, if T is small enough such that C{T,e) ^ i we obtain that 


(25) 

(26) 

(27) 

lL°“(0,T;(S,R^)) oll'^ “ W | |^oo (o,r;(S,K^)) (28) 

This implies that we have a contraction mapping and the hypothesis of 
Banach’s fixed point theorem is satisfied for the mapping 


u 


Il°°(0,T;(S, 


SO that 


^ 2 


1 

^ -||w — w| 1^00(0,r;(S, 


Af : L°°(0,r;(S,M''')) ^ L“(0, T; (S, )) 

W 1-^ 


which therefore has a unique fixed point which solves (fT7|) . If C{T,e) > ^ 
we can choose Ti small enough such that C (Ti, e) ^ ^ and then repeat 
the above argument for intervals [0,Ti], [Ti,2Ti],[nTi,T]. In either case 
we obtain a solution which solves (|17p on the interval [0,T]. Standard 
parabolic regularity theory (see e.g. Th. 5 §7.1 [T7]) then gives us u'^ € 
L^(0, T; 77^(S, R-^)) and that u'" € L^(0, T; L^(S, R-^)), which concludes 
the proof of Theorem [TJ □ 

Note that Evans goes on to use Th. 6 §7.1 [T7] to obtain an improved 
regularity result showing that € L^(0, T; i7^(S, R'^)) and that the time 
derivative € L^(0,T;i7^(S,R^)). However we do not need this result. 
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3.3 Energy estimates 

The next step is to obtain the uniform energy estimate in e for the solutions 
u'". This is the content of the following proposition. 


Proposition 2 

sup (l(^O)-i 


There exists a constant C depending on T, S, hij, dihij, 
, 1^*1, \B\, \dtA^\, \djA'‘\^ such that 


u 


Il2(0,T;L2(S, 


^ C 


|uo||l2(s,r^) + ^ l|f’(^r)llL2(s,RJV)(itj (29) 


Therefore the estimate is independent of e. 

Proof of Proposition\^ Taking the time derivative of I|u^| 1^2(2 gives: 


~ 2 (u^, u'^)^ 2 (s,RAr) 

= 2 (u%eA/,u^ - {A^y^A^diU^ - {A^y^Bu^ + {A^y^p'^ 


(30) 


We have estimates for the following terms in (|3()p 

|(u^,(A°) ^F)l2(s,RJV)| ^ Cl(||W(t)||^2('s 

N 

(W,eAftW)i2(s,RiV) - 
\0\-l 


+ II^IIl2(s 


{u^,iA^y^Bn^) ^ C2||W||i2(s 


= -ej2 I h^^di{y)djiyyh^o 

e \\2 


) (31) 

(32) 

(33) 


However we also require an estimate for (u^, (H^)“^H*9jW)2,2(5] rjv), which 
is the remaining term in (1301) . This term can be estimated by applying a 
suitable integration by parts. 

We hrst assume that the smooth background Riemannian metric is hij = 5ij , 
write {A^y^A^ as and estimate (W, H*5jW)2,2(5] rjv). Using the fact that 
the are symmetric, we then have 


(u*^, HL*9jU’^)^2(s,R:v) (34) 

/s W-(35) 
= ^ ydi •u'^d^x — 5 /s [piA^'^ \T-\TdTx (36) 

= —^ y (^diA^^j u^-u^d^x (37) 

So that 

|(u",i*aiu")i2(s,RA')| ^ ^ ^ W-u"d”x (38) 

^ C'sllu^l 1^2(2,RJV) (39) 
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where the constant is independent of e and we have used the fact that 
diA'‘ is bounded. In the case of a general background Riemannian metric hij 
rather than the ordinary divergence of A*, one obtains the divergence with 
respect to the background metric hij and the corresponding result is 


l(v, 


1 

^ 2 


di + tli)AA 




C4||u^|li2(s,MAf) 


(40) 

(41) 


where again is independent of e. For the case where S is an open, 

bounded set of M” with smooth boundary one needs a slightly more 
complicated argument where one approximates by smooth functions of 
compact support (see m §7.3 for details). 

Using all the available estimates (IHTI) . (I52|) . (1551) and (I4n|) in (IHHIl we obtain 
the estimate 


(||'^^(^)IIl2(s,r^)) ^ ^5 (||u^IIl2(s,r^) + II^IIl2(s,r^)) 

Using Gronwalls’s inequality and the fact that 

l|u^(0,a;)||L2(S,R^) ^ l|uo(a;)||L2(S,R^) (43) 

we obtain the estimate 


sup I|u (i)llL2(S,R^) ^ ^6 
osctsgr ^ ' 



||2 

IIl2(s^rJV) 


+ 



IITTI |2 

11-*^ IIl2(s,RJV) 



(44) 


Finally noting that 


u 


||2 

IIl2(o,T;L2(s,RJV)) 


^ T sup ||u^(t)|| 


2 

L 2 (s,RJV) 


(45) 


and using this in the estimate (|44p we obtain (|29p which concludes the proof. 

□ 


3.4 Existence, Uniqueness and Stability 

In this section we show the existence and uniqueness of the initial value 
problem. In Proposition [1] we obtained solutions u*^ in L^(0, T; i4^(S, M^)) 
to the parabolic system m- Using the Banach-Alaoglu Theorem there 
exists a subsequence that converges weakly to a function u in 

L^(0, T; L^(S,M^)). We now show that this converges to a weak solution of 

(HD- 

First we choose a function w G i?^(S[o_ 7 -],M'^), take the dot product with 
equation (1T7P and integrate over t and x. This gives 
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( 46 ) 


fT 

((A°)■^Lu^w)i2(s,MJV) - e(AfeU^w)i2(s,MJV)(i^ 

~ lo ^^5 w)£,2(s,RAf)C?i 

Integrating by parts we obtain 

(U^, -bw)^2(5]jRAr) — e(u'^, A/jW)j;^2(s,RJV)dt 

~ lo w)^ 2 (s,MAr)dt + (u^(0), w( 0 ))^ 2 (s,RJV) 

where the operator L is dehned by 

Lw := —(9tw — {{A^)~^)-w — f 


(47) 


(48) 


Then taking the limit fc —>■ oo and using the weak convergence of ^ u 
and that u^'“(0) (^'^)“^(0)uo in L^(S,M'^) we obtain 


(u,Lw)i2(s,RJV)dt 


rT 

^ Jo w)j;^ 2 (s,RJV)dt + (uo, w( 0 ))j;^ 2 (s,R^V) 


(49) 


The above equality was obtained for w € however the equa¬ 
tion remains well-defined for w E C C^{[0,T],L'^{T,)) . We 

now show that not only is the equation well-defined, but it remains valid 
for w E -ff^(S[o,r]) J^^)- The method involves taking the convolution with a 
family of mollifiers and then passing to the limit d —>■ 0. 

Let w E -ff^(S[Q 7 i],M'^) and define * w where we have chosen 6 

close enough to zero such that w‘^(T, •) = 0. Therefore we have 



(u, Lw'^)i, 2 (s, 


dt 


(50) 


~ Jo W^)l 2(S,R^)<^^ + (u(0), w‘^(0))i2(s,RjV) 

Taking the limit 5 —)> 0 and using the Schwartz inequality we have the 
following limits 


wf wt in L2(0,r;L2(S,M^)) (51) 

-{A\A^)-^)wf -{A\A^)-^)wim L‘^{0,T-,L^{^,R^)) (52) 

B^{{A°)-^)w^ B^{{A^)-^)w in L^{0,T;L‘^{T.,R^)) (53) 

rff.A''{A°)-^w^ fffc24^(24°)-^w in L2(0,r;L2(S,M^)) (54) 

w'^(O) ^ w(0) in L^(S,M'^) (55) 
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We therefore conclude that 


^ (u,Lw)^2(s,RiV)dt 
rT 

~ ((^'^) w)j;^2(s,RJV)(it + (uq, w(0))j;^2(s,RA') 

for w € 


(56) 


If we now take w G (^““(Sjo^r] j K'^), with supp(w) C [0, T) x S, and multiply 
it by we obtain that G and 

supp(74‘’w) C supp(w) . 

We may therefore insert w = in (|56|1 which gives 


(U, LA°w)^2(s,RA')ht 

— ((^'^) (A°)w)^2(s Rjv)(it + (^uq, (A°)w|t=o) 

which can be rewritten as 


(57) 


L^{T,, 


(u,L*w)i2(s,R^)dt = ^ (F,w)^2(s,RJV)dt+((^°(0))uo,w(0))^^^^^^^ 

’(58) 

for all w G C'°°(Sjq 2 ’],M^) with supp(w) C [0,T) x S, where L* is the 
formal adjoint dehned by equation Q. 

We have therefore proved that the u obtained by taking the limit of the 
subsequence is a weak solution lying in L^(0, T; L^(S, M^)) with 

initial data uq as in Definition [TJ 

As the norm function is lower semi continuous, we may take the limit of 
equation (EHIl to obtain the estimate 

ll^lli2(0,r;L2(S,IR^)) ^ II'^^^IIl2(0,T;I,2(e,RJV)) (59) 

^ C* ^||Uo| 1^2(2, RJV) l|F|li2(s,RA')'^^^ (60) 

So that the solution we have obtained is a regular weak solution. 

To show uniqueness we consider two functions ui, U 2 which are both regular 
weak solutions. Then u = ui — U 2 is a regular weak solution with source 
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F = 0 and initial data uq = 0. Moreover the solution satisfies the energy 
estimate m as shown above. Therefore 

II'^IIl2(0,T;I,2(S,1R^)) = 0 (61) 

This implies u = 0 and therefore ui = U 2 . 

The final step in establishing well-posedness is to prove the stability of the 
solution with respect to initial data. To make the concept precise we say 
that the solution u is continuously stable in L^(0, T; M.^)) with respect 
to initial data uq in L^(S,]R^) , if given e > 0 there is a <5 depending on uq 
such that if uq S with: 


11^0 — ^ 6, (62) 

then the corresponding weak solution u with source function F satisfies 

l|u — u||2,2(o,T;L 2(S,R^)) ^ e (63) 

The stability results follows from using the energy inequality shown in The¬ 
orem \2\ for the difference u — u which gives 

“'^IIl 2(0,T;L2(S,RJV)) ^ C||uo - uoll L2(S,RJV) (64) 

Now choosing <5 = ^ we obtain the inequality: 

11^ “ '^lli2(o,T;L2(s,R^)) ^ (66) 

which establishes stability with respect to the initial data. 

This concludes the proof of Theorem [TJ 


3.5 The wave equation 

In order to apply the results of Theorem 1 to applications in general relativity 
we will show here how the wave equation can be written as a first order linear 
symmetric problem. 

We define 

V = (^lU,..., dnU, dtu, uf = (u\ ..., u”, u’^+2)'^ G M”+2 (66) 

and the symmetric (n -|- 2) x (n + 2) matrices by: 


24° 


^21 „22 


0 0 
V 0 0 


0 0\ 
0 0 

0 

0 l) 
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/ 0 

0 • 

• 9 ^^ 

o\ 

0 

0 • 

• 5^^ 

0 

glk 

g2k . 


0 

\ 0 

0 • 

0 

0/ 


We further define the matrix B to be given by 


0 0 \ 

0 0 

: 0 

y ••• -m2 

V 0 0 ••• -1 0 / 

where are the connection coefficients of the spacetime metric gab- We 
also dehne the M'^-valued vector function F by F = (0,0,... , —/, 0)^. 

In this way, we may rewrite the scalar wave equation m as a first order 
system which has the form 

Lv = A°dtv - A^^^v +Bv = F (67) 

v(0,-) = vo(-) (68) 


/ 0 
0 


B = 


0 

0 


-. ab'pl ^ ab ' p 2 


We may then use the Theorem [T] to establish well-posedness of (I67h . To 
prove existence of solutions to the wave equation, we therefore need to prove 
that the solution v of the symmetric hyperbolic system (I67h has the form 
V = {diu, ..., dnU, dtu, uY. 

We now mollify our solution v using a strict delta net to obtain a sequence 
of smooth functions * v = {v\, ..., that satisfies: 



w)L2(s,RJV)dt 



)L2(E,K^)'^t 


for a suitable mollified F*^ of the form (0,0,... , —/i'", 0)^ where /i'" is any 
function. 

Then 



w)^2(s^KAf)dt — 0 


for all w G C'°°(Sjo,r])l^'^) therefore = F^ almost everywhere so we 
obtain that: 


dtvi = djvY^ 


and 


dtvY^ 


= V 


n+l 
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Also, we choose initial data such that 



( 69 ) 


diu^ir) = •) + 



/ 

Jo 

j 


■T 


■T 


dtv\dt 


didtv^'^'^dt 


(70) 


(71) 

(72) 


Taking now into account that —>■ v in L^([0, T], L^(S, M"'+^)) and that 
the convolution and derivatives commute we obtain the result that v is an 
L^([0, T], L^(S, M”+^)) function of the form {diu, ...,dnU,dtu,u). 

Collecting the results from this section we have established Theorem [2j 


4 Applications 

Although in the following three examples the spacetimes are not spatially 
compact we will assume we are working in a local region of the form S[o,t] = 
[0, T] X T, as described in the geometric setting. 

Junction Conditions There is a precise mathematical formalism proposed 
by Israel to describe the junction conditions for two regular spacetimes joined 
along a non-null singular hypersurface A [7]. He noted that if we consider 
two half-spaces V~^ and V~, a singular hypersurface, A, can be fully char¬ 
acterised by the different extrinsic curvatures (second fundamental forms) 
associated with its embeddings in V~^ and V~ and a continuous matching 
condition of the metric through the common boundary. If we use Gaussian 
coordinates based on A, then the normal derivatives of the metric have a 
jump across A with the metric being continuous along A. This scenario sat- 
ishes the analytic conditions required for the application of Theorem [2] to 
apply. Notice however that the theorem does not need to make assumptions 
on the time dependence or matter content of the spacetime. 

Impulsive Gravitational Waves A spacetime that contains impulsive 
gravitational waves described in double null coordinates has line element 
given by: 


= 2dudv — (1 — uQ{u))‘^dy^ “ (1 + uQ{u))‘^dz^ (73) 
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where Q{u) is the Heaviside step function. The spacetime is vacuum, but 
has a Weyl tensor with delta function components 

Cuyuy — 

Cuzuz — 

It is important to notice that although the curvature is not bounded this 
condition is not relevant for Theorem [2] to apply. 

Brane-world Cosmologies The RandalHSundrum Brane-Worlds (RS) are 
models that explore gravity beyond classical general relativity m. and also 
appear in a cosmological context [52] • In the RS model in AdS^ one has 
that in Gaussian normal coordinates = {t,x^,y) based on the brane at 
y = 0, the model has the line element 

ds'^ = e-2hl/^(-dt2 + dx^^) + dy'^ (74) 

This spacetime again satisfies the conditions for the Theorem [2] to be appli¬ 
cable and therefore solutions with finite energy exist. 

Notice that the well-posedness result can be extended to other brane world 
models and even collision of branes as long as the spacetime satisfies the 
assumptions of Theorem [21 Therefore, one can consider that dynamical 
models of colliding branes (see e.g. [55|) do not produce strong gravitational 
singularities provided that the spacetime remains during all the process. 

4.1 Main result and the relation to previous work 

The main result of the paper (Theorem [2|) shows that spacetimes which are 
usually thought of as singular may be regarded as regular if one adopts the 
point of view that true singularities make the local dynamics of test fields 
ill-defined. We have established general conditions under which linear wave 
equations are locally well-posed in spacetimes with weak singularities where 
the singularity is concentrated on a submanifold. In particular, the results 
can be applied to spacetimes with shell-crossing singularities, surface layers 
and hyper surface singularities of regularity C^’^. 

We establish local well-posedness for general hrst order linear symmetric hy¬ 
perbolic systems with coefficients with low regularity. We show that unique 
stable solutions exist in L^(0, T; L^(S, This solution corresponds in 

the second order formalism to a finite energy solution in of the wave 
equation. Moreover, the main advantage in writing the problem as a first 
order system is that the existence of a covariantly constant timelike vector 
fields and the condition on the curvature are not needed which are key con¬ 
ditions in previous works [Ml IB]- Therefore, the results obtained extend 
previous results of Vickers and Wilson m and Ishibashi and Hosoya M by 
allowing a larger class of non-vacuum time dependant spacetimes. We also 
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establish not only the existence and uniqueness of solutions but also their 
stability and local well posedness. 
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